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Abstract 

According to standard textbooks, the nuclear symmetry energy originates from the 
kinetic energy and the interaction itself. We argue that this view requires certain 
modifications. We ascribe the physical origin of the kinetic term to the discreteness 
of fermionic levels of, in principle arbitrary binary fermionic systems, and relate its 
mean value directly to the average level density. Physically it connects this part 
also to the isoscalar part of the interaction which, at least in self-bound systems like 
atomic nuclei, decides upon the spatial dimensions of the system. For the general 
case of binary fermionic systems possible external confining potentials as well as 
specific boundary conditions will contribute to this part. The reliability of this 
concept is tested using self-consistent Skyrme Hartree-Fock calculations. 
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According to classic nuclear structure textbooks [1] the nuclear symmetry en¬ 
ergy is composed of two basic ingredients: the so called kinetic and interaction 
parts: 


k^sym 2 ^sym 




( 1 ) 


This decomposition and terminology has its roots in the Fermi gas model. 
Within this model the strength of a^n can be easily evaluated to be akin ~ 
100/A MeV. The interaction part, on the other hand, is related to the Hartree 
part of a schematic isospin-isospin interaction: 


Vj'j' — — kT • t. 
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( 2 ) 


The magnitude of aint (or alternatively interaction strength k) can be esti¬ 
mated, however, only indirectly based on the estimate of akin and the actual 
knowledge of the empirical value of the total symmetry energy strength agym- 
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The strength of the nuclear symmetry energy consists of a repulsive volume- 
type XjA part which is reduced at the surface For the purpose 

of this work, however, we adopt an effective volume-like symmetry energy 
strength ttsym ~ 150/^4 MeV. This strength must be understood as the average 
over volume and surface contributions. It is valid for relatively light nuclei of 
yl ~ 50 -h 80 [2], Similar effective values are obtained from hts assuming an 
Esym ~ T(T -|- 1) dependence for the total symmetry energy [3]. Let us point 
out, however, that none of the conclusions drawn below will depend upon the 
adopted strength. 

An alternative way of calculating the symmetry energy is by means of the 
isospin cranking model [4]: 

= H - ujf. (3) 

When the nuclear hamiltonian H can be approximated by an isospin symmet¬ 
ric single-particle (sp) mean-field hgp = '^occ^i fhe model can be strightfor- 
wardly solved analytically. It leads to the isospin parabolic energy law: 

E = (4) 

as shown in Ref. [4,5,6]. The formula (4) was formally derived for an equidis¬ 
tant level model characterized by a sp splitting e. Within such a model the 
sequence of transitions from the T = 0 ground state to T = 2,4, 6, 8, • ■ ■ takes 
place at crossing frequencies which form a simple arithmetic serie e, 35, 5e, • • •. 
In reality, the shell structure will affect that simple pattern leading to a se¬ 
quence of nonuniformly spaced crossing frequencies determined by local shell 
effects. Thus for small values of the isospin strong deviations from the rule (4) 
may show up. However, for large values of T the uniform pattern should be 
recovered with e being the measure of the mean level spacing at the Fermi en¬ 
ergy provided that e{A, Tz) ~ ^{A). Thus, from the perspective of the uniform 
isospin cranking model the ’kinetic’ symmetry-energy term originates rather 
from the mean-level density than the kinetic energy. Note that, apart from the 
smooth mean-part, this term also carries information about the fluctuations, 
directly related to the local shell structure. 


The mean-level spacing in nuclei can be rather easly evaluated. Including 
Kramers and isospin degeneracy it is equal to e = 4:/g{\) where g{\) = 6a/n'^ 
stands for the mean density of states at the Fermi energy. The empirical value 
of the level density parameter a, is uncertain. Typical estimates (again of 
an effective volume-like A-dependent type) range between a ~ A/10 Me 
and A/8MeV“^, see [7,8,9], where the lower limit for a corresponds to the 
harmonic oscillator estimate [1]. Thus the contribution is: 


TT^ TT^ TT^ 53 66 ^ ^ ^ ^ 

e ^ 2— ^ 16—- -=- 20—- ~ —r -i- — MeV. 

3a 3A 3A A A 
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Let us observe that using the Fermi gas model, the estimate for a = A/15 MeV“^ [1] 
in Eq. (5) recovers exactly the kinetic part of the symmetry energy strength 
akin discussed in the introduction. This value of the level density parameter 
is, however, unrealistic small. 

The repulsive isovector term in the nuclear potential (2) can be included within 
the mean-held isospin cranking model. Without loosing generality we can con¬ 
sider the one dimensional isospin cranking model [e.g. around the x-axis, then 
{Ty) = (Tz) = 0]. After linearisation of the interaction (2) the cranking Hamil¬ 
tonian takes the following form: 

^MF ~ ^sp — {ui — k{Tx))Tx ( 6 ) 

with an effective isospin dependent cranking frequency. The isovector potential 
simply shifts the crossing frequencies from e, 3e, 5e, - ■ ■ to e-f k, 3(e-f k), 5(e -t- 
k), ... [within the equidistant level model] leading to 

E = i(E + 2^)T\ ( 7 ) 

The isovector held enters Eq. (7) with a factor of 2fi;, due to the restrictive 
treatment of the linearized interaction as a sp external potential Vt = k{T)T. 

In this manner, the isovector term is taken into account in phenomenological 
potentials. Although n can be adjusted to yield the proper Fermi energy de¬ 
pendence on T, summing up the total energy from the occupied sp energies, 
results in a gross overestimate of the symmetry energy. Hence, estimates of k 
based on phenomenological potentials are not very reliable. Indeed, only by 
means of the Strutinsky shell-correction procedure one is able to calculate the 
total energy. This implies, however, that the information carried out from the 
potential is stripped oh any average trends including the mean level density 
and K in particular and only the huctuations are transfered. 

Nevertheless, the Strutinsky smooth energy can be used for crude verihcation 
of the reliability and stability of our concept. We calculated the smooth Struti¬ 
nsky energies for a chain of e-e A=88, = 0 -i- 10 nuclei. We increased the 

strength of the isovector potential in steps from zero to its full value. It turned 
out, that the level density, obtained in the calculations does not depend on 
Tz- This suggest that this method can be viewed as an independent, accurate, 
and simple way of calculating the mean level density at the Fermi surface. 

Within the Hartree-Fock (HF) approximation, the cranking model (Eq. (3)) 
including an interaction Vtt (Eq- (2)) will result in: 

E =-{e + k)T‘^ +-kT . (8) 

Thus, the interaction (2) leads to a linear term which comes entirely from 
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the Fock exchange [10] or, in other words, results in Egym ~ T(T + x) with a 
X < 1 dependence for the symmetry energy. This form of the symmetry energy 
is intermediate between Egym and Egym ~ T(T + 1) like the formulas 

debated in the literature. Interestingly, the self-consistent HF or Hartree-Fock- 
Bogolyubov (HFB) models are believed to give a quadratic Egym ~ type 
dependence although, to the best of our knowledge, this conjucture has never 
been discussed seriously. A dependence ~ T(T -|- 1), on the other hand, can 
be deduced from the nuclear shell-model. 


The interaction (2) is highly schematic as compared to the isovector part 
of e.g. a Skyrme type interaction.^ Therefore, the symmetry energy formula 
(8) can serve only as a guideline for further investigations which are entirely 
based on the Skyrme-Hartree-Fock (SHF) energy functional. Nevertheless, we 
will treat it literally and show that the entire concept of strictly dividing the 
symmetry energy into contributions from mean-level-density and interaction 
is very reliable and highly transparent also from the point of self-consistent 
helds. To investigate the isovector effects of the strong force we switched off 
the Coulomb interaction and assumed the equality of proton and neutron 
masses. The calculations presented below were done using SHF code HFODD 
of Ref. [11] employing various parametrisations of the Skyrme forces. 

The mean-lev el-density e contribution can be studied simply by removing the 
isovector part of the Skyrme force. Then, guided by Eq. (8), the calculated 
HF energy is approximated by: 

AB*""' = Bt""’ - (9) 

which allows to extract e as a function of T(= T^). The results of the calcu¬ 
lations are presented in Fig. 1. The calculations have been done for several 
Skyrme forces including SLy4 [12], SHI [13] and SkO [14] i.e. forces having 
isoscalar effective masses equal m*/m ^ 0.70,0.76 and 0.90, respectively, as 
well as SkP [15], SkXc [16] and MSk3 [17] i.e. forces having m*/m ^ 1. The 
calculations were done for the even-even A = 48, 68 and 88 isobaric chains of 
nuclei from = 0 to the vicinity of the drip line. 

As anticipated, for relatively small values of strong variations in e are seen. 
However for larger values of > 8 the values of £ stabilize. Let us point 
out that the forces having effective masses considerably lower than unity have 
a substantially larger £. However, after employing an effective mass scaling 
e —>■ all curves are within the limits given by the estimate (5) [shaded 

^ To be precise, in the following, the isovector part of the Skyrme interaction always 
refers to that part of the interaction that gives rise to the isovector dependence of 
the mean nuclear potential. Within the standard Skyrme energy density functional 
it is given by the isovector density dependent terms. 
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Fig. 1. The mean-level spacing (left hand side) calculated using only the isoscalar 
part of the SHF interaction for the A = 48, 68, 88 isobaric chains. Parametrisations 
used are indicated in the legend. Right hand side shows the results scaled by the 
isoscalar effective mass Shaded areas indicate the estimates given by (5). 

areas in the figure]. Clearly, the physical interpretation of this term as arising 
from the level density is evident. This observation nicely conhrms the validity 
of our approach. Let us further observe that, with increasing A, all curves 
move towards the upper limit of (5). This, most likely, reflects the decreasing 
role of surface effects with increasing A. Indeed, more accurate estimates of 
the average density of states g, including surface effects can be done based 
on the semiclassical theory [18]. It can be shown that surface effects increase 
the density of states gv+s > 9v and, in turn, lower average level spacing 
Ey+s < thus giving rise to the observed trend. A detailed discussion of the 
A dependence of the symmetry energy based on these ideas will be published 
elsewhere. 

The dependence of e on the effective mass is directly related to the isovector 
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Fig. 2. The average effective strength k of the isovector part of various 
Skyrme forces. Left hand side shows the values of k estimated assuming 
Esym = (^ + E)T’^I2. Right hand side shows the values of k estimated assuming 
Esym = !2 + kT(T + 1)/2. In both cases the values of e from the isoscalar SHF 
calculations have been subtracted. 

part of the Skyrme interaction. A small effective mass results in a reduction 
of the isovector Skyrme interaction, since the overall symmetry energy must 
reproduce, by construction, the empirical value. This effect is seen in Fig. 2 
showing the estimated values of k which, in this context, have the meaning 
of an effective strength of the isovector component of the Skyrme force. We 
have performed two set of calculation to extract k. In the calculations pre¬ 
sented in the left panel of Fig. 2, we have assumed a purely quadratic isospin 
dependence = |(e -f k)T‘^. The right panel, on the other hand, shows 

calculations that include also the linear term from Eq. (8). For both cases we 
have made use of the e values extracted previously (Fig. 1) from the isoscalar 
SHF calculations.The calculations reveal at least two striking features: {i) No 
major perturbations in the k{Tz) curves are seen for large pointing to a 
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very weak, if any, impact of the isovector fields on the average level density, 
(u) Clearly, and rather surprisingly, these calculations reveal the presence of 
a linear term. Indeed, the curves in the right panel of Fig. 2 show essentially 
no Tz dependence. The only exception from this rule is the SkO force with its 
highly non-standard strong, repulsive isovector component of the spin-orbit 
force. 

An interesting conjecture can be made here concerning recent very accurate 
self-consistent mass calculations based on Skyrme forces. [19] These calcu¬ 
lations clearly prefere Skyrme forces of m*/m ~ 1. This preference can be 
(partly) related to the presence of the linear term in the symmetry energy, 
provided, that the empirical symmetry energy has a term like ~ T(T -|- x). 
Indeed, for Skyrme forces having m*/m < 1 the linear term is expected to 
be reduced. On the other hand, there seems to be no counterpart to the 
~ |eT type at the mean-held level. Such a compensating term may be calcu¬ 
lated within the RPA theory as suggested recently by [10]. Performing RPA 
calculations atop of HF/HFB particularly in large-scale mass calculations is 
presently beyond reach. However, since the linear term describes, in fact, the 
isospin dispersion, and since our interpretation links the inertia parameter to 
the average spacing e, the overall correction can easily be estimated based 

on Eq. (5) [scaled by the isoscalar effective mass] and the calculated isospin 
^2 

dispersion (AT ). 

In the sp limit, one hnds that (AT^) = T and therefore, the interaction part of 
the symmetry energy follows ~ T(T -|- 1). However, pairing correlations 
tend to increase (AT^), with the largest effect in N=Z nuclei, as shown in 
the insert of Fig. 3. Thus pairing effectively weakens the linear term leading 
to rsj T(T -|- x) with a; < 1. The impact of pairing on the value of x is 
illustrated in Fig. 3, and dehned as 

^^(ATV^^AT^^ where T = T, = (A - Z)/2 (10) 

The value of x is calculated using the schematic, non-selfconsistent constant 
gap A BCS calculations. In the calculations the equidistant level spectrum 
characterized by a sp splitting e was used. 2N (2Z) of the double degener¬ 
ate levels were active in the pairing calculations, in order to take advantage 
of the particle-hole symmetry which in turn guarantees automatic number 
conservation. 

As already discussed pairing has the greatest impact on A ~ Z nuclei reducing 
the value of a; to ~ 0.3 for the A = 88, = 1 system, see Fig. 3a. With 

inreasing T, the value of x(T) increases slowly to a; —1. For small T the 
actual value of x changes rapidly with an increasing ratio of A/e as illustrated 
in Fig. 3b. Since A/e ~ [z/ ~ 1/2 -h 2/3] the quenching of x is expected to 

be more efficient in heavy N Z nuclei. Local shell-structure, however, may 
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Fig. 3. BCS calculations, based on a schematic, equidistant level model of the pa¬ 
rameter X (Eq. 10) in the interaction related part of ~ kT(T -|-x). The upper 
part shows x{Tz) calculated for A = 88, A = 12/\/)dMeV and two values of e as 
indicated in the figure. The lower part shows the quenching of x versus A/e for 
Tz = 1 [the filling approximation is used], = 2 and = 10 nuclei. The mecha¬ 
nism responsible for the quenching of x in A ~ Z nuclei is displayed in the inset 
which shows the T-dependence of the interaction related linear term k(AT^)/ 2 for 
A =0 (single-particle) 1, and 2MeV, respectively. 

result in oscillations of the A/e ratio and, in turn, to strong variations in x. 

To verify this assumption we performed self-consistent HF-I-BCS calculations 
using SLy4, Sill, SkP and MSk3 forces. The calculation scheme was similar to 
the one described above. First we extracted Sbcs from the isoscalar Skyrme 
SHF-I-BCS calculations. Then, using these values, we extracted kbcs- The 
values of kbcs were computed assuming ~ kT{T + 1) and compared 

to corresponding values from pure SHF calculations, khF) see Fig. 2. The 
difference is shown in Fig. 4. It is always negative, and particularly large for 
small Tz- This clearly conhrms the schematic model calculations shown in 
Fig. 3. 
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Fig. 4. The difference between HF+BCS and HF values of k. In both cases k has 
been calculated assuming ~ kT{T + 1). Note that values of kbcs — i^hf 

are always negative. This behavior shows that pairing has a similar destructive 
impact on the linear Skyrme related term as emerged from the simple schematic 
considerations shown in Fig. 3. 

In conclusion, it is demonstrated that the Skyrme forces give rise to a linear 
term in the interaction part of Egym which behaves effectively similar to the 
schematic interaction (2), i.e. like ~ (AT^). This term is strongly quenched 
in TV ~ Z nuclei due to pp/nn pairing, which, even after including RPA 
correlations [10], seem to leave room for isoscalar pairing correlations mN = Z 
nuclei. Particularly, since the data clearly indicates that the linear term of the 
symmetry energy in the vicinity of the N = Z line is enhanced and goes 
like ~ T(T + 1.25) [20,21]. Mass calculations with the SHF need to take into 
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account the fluctuations in isospin, ~ (AT^) similar to the fluctuations in 
spin [17] with an inertia given by the mean level spacing at the Fermi energy. 
This may allow for more general Skyrme forces, not restricted to effective mass 
m* ~ m. 

In particular, the present investigation reveals that one component of the nu¬ 
clear symmetry energy can directly be associated with the mean-level density 
rather then the kinetic energy as proposed in textbooks. Although the argu¬ 
ments are based on Skyrme-Hartree-Fock calculations of A = 48, 68, and 88 
chains of nuclei, this is a general feature originating from the granularity of 
the fermionic single particle spectra and not related to atomic nuclei only. 
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